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^ ■ Abstract 

In this article we present calculations which indicate that nanoshells can be used as a high 
pressure gauge in Diamond Anvil Cells (DACs). Nanoparticles have important advantages in 
\ comparison with the currently used ruby fluorescence gauge. Because of their small dimensions 

they can be spread uniformly over a diamond surface without bridging between the two diamond 



^ ' anvils. Furthermore their properties are measured by broad band optical transmission spectroscopy 

leading to a very large signal-to-noise ratio even in the multi-megabar pressure regime where ruby 
measurements become challenging. Finally their resonant frequencies can be tuned to lie in a 
j3 \ convenient part of the visible spectrum accessible to CCD detectors. Theoretical calculations for 

a nanoshell with a Si02 core and a golden shell, using both the hybridization model and Mie 
' theory, are presented here. The calculations for the nanoshell in vacuum predict that nanoshells 

^ . can indeed have a measurable pressure-dependent optical response desirable for gauges. However 

O . when the nanoshells are placed in commonly used DAC pressure media, resonance peak positions 

as a function of pressure are no longer single- valued and depend on the pressure media, rendering 
them impractical as a pressure gauge. To overcome these problems an alternative nanoparticle is 
' studied: coating the nanoshell with an extra dielectric layer (Si02) provides an easy way to shield 

. the pressure gauge from the influence of the medium, leaving the compression of the particle due 

, to the pressure as the main effect on the spectrum. We have analyzed the response to pressure 

\ up to 200 GPa. We conclude that a coated nanoshell could provide a new gauge for high-pressure 

measurements that has advantages over current methods. 
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I. INTRODUCTION 



Recently nanoshellsjH have received a lot of attention because they possess interesting 
optical properties. Consisting of a dielectric core surrounded by a metallic shell, nanoshells 
allow surface plasmon polaritons to exist at the interfaces between the layers, resulting in a 
spectrum with very pronounced absorption and extinction peaks. The position and broad- 
ness of these resonance peaks depend on the size, shape and materials of the particles and on 
the surrounding medium. This versatility and sensitivity make nanoshells excellent probes in 
different fields, ranging from cancer ablation [2| and biosensors 0] to surface enhanced raman 
spectroscopy applications j3| and plasmonicsjH]. In this paper we investigate the usefulness 
of nanoshells in the field of high-pressure physics, and more precisely Diamond Anvil Cell 
(DAC) pressure measurements [6(]. The first theoretical calculations for a spherically sym- 
metric Si02-gold nanoshell under pressure (up to 200 GPa) are presented, indicating clear 
measurable absorption peaks. The main idea behind the use of nanoshells as a pressure 
gauge is that the size of the nanoshell and the permittivity of the Si02 core and medium 
will change under pressure resulting in a shift of the resonance frequencies. The goal of 
this paper is to calculate the optical transmission as a function of wavelength for a spatial 
distribution of noninteracting monodisperse nanoshells. 

At the moment there are several pressure gauges used in high-pressure physics. We 
briefly discuss them and conclude with the advantages of the nanoshell-gauge in comparison 
with the established methods. The prime and venerable pressure gauge is the shift of the 
peak of the ruby fluorescence spectrum with pressure 0, Hj, excited by laser pumping. In 
this case micron sized ruby grains are embedded in the pressurization medium. Several 
problems can arise above a certain pressure jo], lo[: at pressures above 100 GPa diamond 
fluorescence can become intense and mask the ruby line; ruby grains can be blown away 
during loading or masked during pressurization; above a few hundred GPa the standard 
method of pumping the ruby with a green or blue laser line becomes ineffectual as the 
laser light is absorbed by the diamond anvil. Although some of these problems can be 



overcome I 111. |12| . ruby becomes challenging to use at very high pressure. For cases where 
the ruby is lost or cannot be observed researchers have used the phonon Raman spectrum 
of the diamond from the high pressure culet or stressed region of the diamond. This is 
less precise and has some dependence on the diamond geometry (l3|. Another method is to 
embed a grain of diamond in the pressurization medium, but this also has challenges and 
limitations |l1^]. An important gauge is the X-ray spectrum of metal "markers" embedded 
in the pressurization medium, but this is only useful at synchrotrons. Finally we mention 
that many of the pressure gauges have problems at high temperatures. The advantages of 
the nanoshells, which come dispersed in a volatile fluid, is that they can be painted onto 
the diamond culet in a thin invisible (to the eye) layer so they will not blow out and they 
cover the entire culet flat so they cannot easily be masked; they should maintain their 
sensitivity to the highest pressures; optical spectroscopy is easy to implement with a large 
signal-to-noise ratio, compared to fluorescence or Raman scattering, and will not be masked 
by fluorescence from the diamonds. We believe that nanoshells will maintain their sensitivity 
at high temperatures, but this requires further study. 

The article is organized as follows. In section II we first review existing theories that 
can be used to calculate the spectrum and resonance frequencies of a nanoshell. Next, the 
pressure dependency is taken into account. In section III the results of the calculations are 
shown for the nanoshell geometry. To overcome some of the problems seen with nanoshells 



2 



the coated nanoshell is introduced in section IV. Section V contains the conclusions. 



II. THEORY 

The optical properties of nanoshells are usually described by either Mie theory or the 



hybridization model. The hybridization model, developed by Nordlander and co-workers |1 5 
offers a fast-to-evaluate analytical form for the resonance frequencies, but it is only valid for 
small particles (smaller than ca. a tenth of the wavelength). Mie theory [l6| is valid for all 
particle sizes but this method is far more time consuming for calculating the peak positions 
since it does not allow for a closed form analytical formula for the position of the resonance 
peak. These resonance peak positions have to be numerically derived from the absorption 
or extinction spectrum. 

In this section first both theories will be reviewed. The pressure dependency is taken into 
account in the last part of this section. 



A. Hybridization theory 

Hybridization theory[15'] is a phenomenological approach in which the nanoshell is mod- 
eled as a combination of a cavity and a metallic sphere. First the sphere and cavity are 
studied separately. The free electrons of the sphere are descibed as an incompressible fluid 
excited by plasmons and the dielectric function is assumed to be that of a Drude metal[l7|: 

Es (w) = £6 - (1) 

CO 

where Eb indicates the background permittivity and cOpi the bulk plasma frequency of the 
metal. 

In a second step the resonance frequencies of the cavity and the sphere are hybridized due 
to the coupling between the cavity and the sphere modes, very similar to the hybridization 
of molecular orbitals. This results in a closed formula for the resonance frequencies jiol]: 

= ^{{L + if Em + 2L (L + 1) Eh + Lh, + L (L + 1) {e^ - 2^^ + e,) x^^^^ (2) 
± { [Lh, - (L + 1)' + L (L + 1) {Em - £c) x^"^^^] ' 
+4L (L + 1) [{L + l)Em + LE^f x^^+i}'/') 

with: 

/3 = 2 { [(L + 1) + Lec] [{L + l)Em + Lei] (3) 

+L {L + 1) {Eb - E,) {Em - Eb) x'^+l} , 

where Ec, Eb and Em are the permittivities of the core, the background of the shell (see 
formula ([1])) and the medium respectively, L is the orbital quantum number and x = Rc/Rs 
is the ratio between the core (i?c) and shell {Rg) radii. 
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TABLE I: The material properties used in the pressure calculations of the nanoshell. Parameters 
Kq and Ki are valid up to 580 GPa for gold and 150 GPa for amorphous silicon according to resp. 
ref. a and ref. 0. The last column indicates the relative permittivity at ambient pressure. 



The permittivity indicated for gold is the bulk permittivity as defined in equation [TJ 



Material 


Ko (GPa) 


Ki 


e 


Amorphous Si02 
Au 


329^22] 
167f24] 


4.1 [22] 

5.94[24] 


2.0449i22J 

6.9[25J 



B. Mie theory 



Mie theory is used to solve the Maxwell equations in systems with spherical symmetry [16 
It allows the calculation of the absorbed, scattered and total cross section for an incident 
beam of radiation as a function of frequency. The nanoshell is modeled as two concentric 
spherical spheres. For this case the results originally derived by Mie for a single spherical 
interface were worked out by Aden and Kerker(20|. The coated nanoshell, which will be 
introduced later, is modeled by three concentric spheres. Extending Mie theory to three 
concentric spheres, although not available in literature, does not provide any new problems 
in comparison with the two concentric sphere case solved by Aden and Kerker. The final 
expressions for the electromagnetic fields in and around the coated nanoshell, as well as the 
scattering and absorption cross section, are given in the appendix. From the resulting total 
cross section the resonance frequencies can be determined numerically. 

For this method the input parameters are again the sizes and permittivities of the shells. 
In analogy with the hybridization model we used the Drude dielectric function in Mie theory 
to describe the permittivity of the metal shell. 



C. Pressure dependency: equation of state 



Applying a pressure P to a nanoshell will lead to a change of the volume of the core 
Vc (P) and the shell Vg (P). These changes can be calculated using the Vinet equation of 
state (E0S)|2l|: 



P = 3K, 



' Vo 



0" 



1/3 



Vo J 



2/3 



exp <j - (iTi - 1) 



1/3- 



(4) 



which gives the relation between the pressure P and the relative volume change V/Vq with 
respect to a reference volume Vq of a material at a certain temperature. Here Kq and Ki are 
material constants for which the values used in this article are presented in table [H Using 
the Vinet EOS one can calculate the size change of the nanoshell under pressure where Vq 
has been taken as the original volume at zero pressure. Under pressure the nanoshell will 
be compressed resulting in a different core and shell radius and a different bulk plasma 
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FIG. 1: The behavior of the sizes of the nanoshell under pressure. The top panel presents the radii 
of the core and shell as a function of pressure. Notice that the golden shell will be compressed more 
than the core until about 150 GPa where the core/shell ratio reaches a maximum. The bottom 
panel shows this ratio Rc/Rs, which is an important parameter for the optical response. The radii 
of the nanoshells in this article correspond to the size of commercially available nanoshells. 



frequency. These changes can be calculated from: 



K(0) 



K(o) 



K(0) 



K(0) 



1/3 



n (P) 



VjP) ) ""-^ 



(5) 

(6) 
(7) 



where Vc (P) and K (P) are the volumes of the core and the shell respectively, n is the 
bulk electron concentration, e is the electron charge, Eq is the permittivity of vacuum and 
m is the electron mass. The pressure dependency of the radii is shown in Figure [1] for a 
nanoshell with = 60 nm and Rg = 75 nm. When the nanoshell is compressed the electron 
density will increase and so will the bulk plasma frequency dTj). This will also influence 
the permittivity of the metal since the plasma frequency is an essential part of the Drude 
dielectric function (see Eq. ([1])). 

The pressure and frequency dependency of the permittivity of Si02 in the megabar regime 
remains, to the best of our knowledge, unknown. However it is possible to estimate the 
pressure dependency from the Vinet EOS using the Clausius-Mossotti relation |l7l|: 



ei^,P)-l 471 

= — n (P) a, 

e{u},P) + 2 3 ^ ^ ' 



(8) 



where n (P) = N/V (P) is the electron density and a the average atomic polarizability, 
assumed independent of pressure. This relation links the permittivity to the volume of 
a material. Dividing this equation with the same equation at zero pressure allows us to 
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FIG. 2: The pressure dependent permittivity of different materials. Tlie permittivity of Si02 was 
calculated from the Vinet EOS and the Clausius-Mossotti relation. The other permittivites are 
based on literature 27f| . 



estimate 

VciO) e(w,0)+2 

where V (P) /Vq can be calculated from the Vinet EOS (jlj) and e {u, 0) represents the relative 
permittivity at zero pressure. The pressure dependency of the permittivity of Si02 and some 
other possible medium materials is shown in Figure |2l Notice that there is no phase change 
in the Si02 curve. In nature Si02 is in a crystalline phase at zero pressure and will switch 
to the amorphous phase at about 90 GPa. However, it has been reported [26[ that the 
fabrication method of the Si02 particles results in cores which are in the amorphous state 
also at zero pressure. 



III. RESULTS FOR A NANOSHELL 



Figure |3] shows the relative cross section of the considered nanoshell placed in vacuum 
as calculated with Mie theory. The relative cross section is the total optical cross section 
divided by ttR'^, the area of the projection of the nanoshell on a plane perpendicular to the 
incoming radiation. From this figure one can clearly see a substantial blueshift when 200 GPa 
pressure is applied. Another important observation concerns the width of the peak. This is 
important for increasing the precision in an experimental measurement since sharper peaks 
allow for a more accurate determination of the peak position. However, sharp resonance 
lines may elude experimental detection if they do not carry enough spectral weight. In 
Figure [3] one can see that the dipole peak (the rightmost peak) is a broad, clear and rather 
symmetric peak, while the quadrupole peak is much sharper and a good candidate for these 
kind of experiments. Numerically however it is easier to track the broad dipole peak, thus 
all results presented here will be with regard to the dipole peak. It is seen that although 
the quantitative results differ for all peaks, the qualitative results presented here hold true 
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FIG. 3: The relative cross section as a function of wavelength for different pressures for a nanoshell 
with parameters as mentioned in the figure. Notice the pressure-induced blueshift of the dipole 
peak from 611 nm at GPa to 531 nm at 200 GPa as indicated by the two black lines. 
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FIG. 4: The position of the dipole resonance peak as a function of the pressure for a nanoshell 
with parameters as mentioned on the figure. The two horizontal black lines correspond to the two 
black lines on figure El The blueshift from 611 nm at GPa to 531 nm at 200 GPa as predicted 
by Mie theory is clearly visible. The red circles were calculated with hybridization theory and also 
indicate a blueshift, but predict a different peak position. The inset shows the peak shift as a 
function of pressure compared to the original position at zero pressure. 



for all resonance peaks. 

In Figure jH the position of the dipole resonance peak is shown as a function of pressure 
for a nanoshell in a medium with e = 1 constant. The black squares are calculated using Mie 
theory and a subsequent numerical determination of the peak maximum. The results from 
the hybridization theory are presented by the red circles. It is clear that the two theories 
do not agree on the exact position of the resonance peak for the considered nanoshell. This 
is due to the electrostatic limit used in the hybridization model which is only acceptable for 
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FIG. 5: The position of the dipole peak as calculated with Mie theory for various pressure media: 
helium [2 7l| (red circles), hydrogen [28| (blue diamonds) and neonj27|| (green triangles). In all figures 



the vacuum position is indicated in black for reference and the phase transition from liquid to 
solid is indicated. The plots only show the pressures for which experimental data is available. For 
completeness the original data are presented in table HI] together with the pressure range in which 
they are valid. 



nanoparticles much smaller than the wavelengths of the incident light (the limiting case for 
nanoshells much smaller than the wavelengths used in a DAC experiment, with sizes at the 
moment unachievable for fabrication, do agree). Still the hybridization model can prove to 
be useful because both theories agree well on the amount the peak shifts, as is shown in 
the inset of Figure |H The two theories diverge from each other only at higher pressures. 
The results show that the nanoshell could be used to measure the pressure by determing 
the amount the resonance peak shifts. The inset is also an indication of the resolution that 
could be achieved. For this example nanoshell the dipole resonancepeak shifts over 75 nm 
when 200 GPa pressure is applied (for rubies this would be 51.1 nmj^, [lo|). 

To determine the usefulness of nanoshells as a pressure gauge, calibration would need to 
be done for every pressure medium separately. The nanoshell's optical response depends 
on the dielectric function of the surrounding medium in which it is placed, and as such 
both the position of the resonance peaks and their pressure-induced shift will differ for 
different pressure media. The original data used for the calculation of the dielectric function 
are shown in table [Tll In Figure [5] one can see the pressure dependency of the position of 
the dipole resonance peak for different pressure media: helium, hydrogen and neon. The 
behavior of the medium under pressure clearly has a large influence on the optical response 
of the nanoshell. For these media, the blueshift that occurs in vacuum has turned into 
a redshift as pressure is increased. From a certain pressure on, there will be almost no 
shift of the resonance peak, rendering the nanoshells ineffectual as pressure gauges in that 
region. Although it is possible to calculate the optical response of nanoshells for much higher 
pressures, the results are not presented here because they are based on extrapolation of data 
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TABLE II: The refractive index of materials used as pressure medium in a DAC as reported by 
references 27| and 28|. The last column presents the pressure range of the data on which these 



fits are based. The dielectric function e can be calculated by squaring the refractive index. 



Material 


Refractive index 


P range (GPa) 


He (fluid) ^27] 


n = 0.8034 + 0.20256 (1 + p)0-i2763 


0.08 - 11.5 


He (solid) aa 


n = -0.1033 + (1 + P)°°^2 


11.7-20.2 


H2i28J 


n = -0.687343 + 0.00407826P 
+1.86605 (0.29605 + p)0 0646222 


0-100 


Ne (fluid) ^27] 


n = 0.668 + 0.33 (1 + 4.3P)°-°'^^ 


0.7-4.7 


Ne (solid) f27] 


n = 0.9860 + 0.08597P°-i953 


5-27 



from a limited pressure region and as such are considered unreliable. 



IV. RESULTS FOR A COATED NANOSHELL 

As was seen in the previous section, the simple nanoshell geometry is not ideal for high- 
pressure experiments. The peak shift is affected by the pressure medium so calibration will 
be necessary for each nanoshell and medium. A possible solution would be to shield the 
nanoparticle from the effects of the medium, therefore allowing the shift of the resonance 
peak to be influenced only by the compression of the nanoparticle. The absence of these 
problems when the nanoshell is placed in vacuum indicates that shielding the particle could 
indeed result in the desired effect. 

Shielding the nanoparticle from the environments means creating a barrier between the 
golden shell and the dielectric medium. This barrier should also be a dielectric. If not, an- 
other set of surface plasmons polaritons would arise on the interface between the outer layer 
and the environment, counteracting the intent of the extra layer. In this article the extra 
coating is achieved by adding an extra Si02 layer to the model and extending Mie theory 
to three concentric spheres. Henceforth we shall indicate this type of " nanomatryushka" [iSj 
as "coated nanoshell". The main question to be answered is how thick this coating should 
be to effectively shield the nanoshell from the pressure effects on the dielectric function of 
the environment. Figure E] presents the dipole peak position for several coating thicknesses 
of the coated nanoshell as calculated by Mie theory. For reasons of comparison the core and 
the thickness of the golden shell are kept the same as for the nanoshell discussed before. 
Furthermore we have only studied the coated nanoshell in a helium environment since this is 
a commonly used quasi-hydrostatic pressure medium in DAC experiments. In neon the re- 
sults will be almost the same because the dielectric functions of helium and neon are similar. 
For hydrogen the counteracting effect is not present, therefore the results with nanoshells 
are adequate and no further improvements are necessary. 

In this article we do not compare the calculations for the coated nanoshell with calcula- 
tions using the hybridization model. The expansion of the hybridization theory for multiple 
layered nanoshells is available in literature [l9|, but it is not considered here due to the ex- 
pected deviations since the diameters of the coated nanoshell particles are 2 or 3 times larger 
than the nanoshell particles and thus the electrostatic approximation is certainly not valid. 

It can be seen from Figure [6] that for thin coatings the redshift due to the effect of the 
medium is still visible. At a certain pressure the peak position will reach a maximum and 
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FIG. 6: The dipole peak position as a function of pressure for the coated nanoshell with a 60 nm 
Si02 core, a 15 nm thick golden shell and for different thicknesses of the outer Si02 coating. For 
small coatings the redshift due to the medium is still clearly visible, while for thick coatings this 
effect seems to disappear. The black crosses indicate the pressure at which the wavelength peak 
position is maximal and where the redshift turns into a blueshift. 

from there on the peak will undergo a blueshift which is due to the compression of the 
nanoshell similar to the results in vacuum (Figure H]). This maximum is indicated by the 
black crosses and can be used as a measure for the pressure up to which the medium affects 
the optical response of the nanoparticle. It is clear that for thicker coatings the maximum 
shifts to lower pressures suggesting that the influence of the medium indeed diminishes and 
thus that the coating effectively shields the nanoparticle. Unfortunately the redshift never 
disappears completely meaning that the influence of the medium cannot be fully shielded by 
the dielectric layer. However, it is possible to position the maximum into a pressure region 
where it can do no harm for the pressure measurements. 

Obviously, the thickness of the coating has to be adjusted to the needs of the experiments. 
A thicker coating will provide better shielding, but bigger particles take up more space and 
risk interacting with each other due to being in close proximity of one another. However, 
since the coated nanoshell is much larger than the nanoshell, the absolute cross section is 
larger, meaning that a smaller density is sufficient while still gaining the same response. 

The derivative of the curves on Figure E] give a direct indication of the sensitivity achiev- 
able with the coated nanoshell structures. On average the sensitivity for the coated nanoshell 
with a 110 nm thick coating (upper curve on Figure [6]) is approximately 0.90 nm/GPa. 
The resolution for rubies in the same pressure regime (0 to 20 GPa) is approximately 
0.36 nm/GPa^ and for a nanoshell in vacuum as was shown in Figure H] this would be 
0.94 nm /GPa. From this we can conclude that the theory predicts a better resolution for 
nanoshells than for rubies in the pressure regime under consideration. 

V. DISCUSSION AND CONCLUSIONS 

In this article we have carried out a theoretical analysis of nanoshells which can be 
designed so that they have absorption peaks in the IR-visible part of the optical spectrum 
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due to scattering or absorption by localized surface plasmon polaritons. Calculations indicate 
that for nanosliells in vacuum these peaks have a substantial shift with pressure making them 
suitable as a pressure gauge for high-pressure research. A useful pressure gauge should have 
a calibration (peak wavelength vs pressure) independent of the pressurization medium. We 
found that for a simple nanoshell consisting of an Si02 core and a gold shell, the calibration 
differed with the pressurization medium; moreover it was double-valued and had a region of 
zero slope. The latter problem was resolved by coating the nanoshell with an Si02 cladding, 
resulting in a robust sensitive pressure gauge. The proposed nanoshell gauge has advantages 
in comparison with well-established pressure gauges: they will not easily blow out during 
loading of a DAC, they cannot easily be masked by fluorescence from the diamonds, they 
have large signal-to-noise ratio in comparison to fluorescence and Raman Scattering, and 
due to their small dimensions they will not be stressed by bridging between the diamonds 
at very high pressure when the gasket thins. All of these putative advantages must be 
confirmed by experiment. 

A possible implementation of the coated nanoshell would be to distribute them inside 
the DAC cell together with the ruby on the diamond culet. For low pressures both gauges 
can be used and the coated nanoshell can be calibrated by using the extensive knowledge of 
the behavior of ruby under pressure Q]. A possible challenge is the application of nanoshells 
to the surface of a diamond culet. Coated nanoshells can be acquired at high concentration 
(2 ~ 3 X 10^/ mm^) from Nanospectra Biosciences, Inc., in a liquid solution. A droplet can 
be placed on the culet and allowed to evaporate to produce a coverage bonded to the surface 
by van der Waals forces. Preliminary measurements show that to avoid clustering and 
segregation, it may prove useful to functionalize the diamond surface with a film of poly- 
4-vinylpyridine (pvp) which has dense sites that localize the nanoparticlesfsO^. For high 
pressures the ruby measurement would be difficult or no longer be possible. The spectra of 
the coated nanoshells however will still be measurable since these measurements are based 
on absorption and transmission. In this way nanoshells could be easily used and effectively 
extend pressure measurements to ultra high pressures. 
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Appendix A: Absorption and scattering of an electromagnetic wave by 3 concentric 
spheres. 

The coated nanoshell was modeled using Mie theory for 3 concentric spheres. These 
spheres divide the space into 4 regions: the core of the nanoparticle, two shells and the 
medium in which the nanoparticle is immersed, respectively numbered 1 through 4 as shown 
in Figure [71 The radii of the three spheres are indicated with a, b and c, from inner to outer. 
Using the notations of Bohren and Huffman [29| the incoming plane wave can be expanded 
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Region 4 




FIG. 7: Schematic overview of the modeled system: region 1 will be the Si02 core, region 2 the 
golden shell, region 3 the Si02 coating and region 4 the pressure medium. 



in vector spherical harmonics: 



oo 

E {^oll - ^N;,ll) , (Al) 



n=l 

, oo 
— Kd \ — > „ / . 



CU/i4 

n=l 



H.,^ = :^J2^AMSn + ^N:.ll), (A2) 

with: 



En = ^0^"4^- (^3) 

n{n + 1) 



The vector spherical harmonics can be written as: 
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On the right hand side one can see the wel known spherical Bessel, Neumann and Hankel 
functions respectively. 

In region 4 there will also be a scattered wave which can be described by: 

OD 

E,,s = Y.^n {-anMSi + ^KN^li) , (A13) 

n=l 
, oo 

H4,s = ^ ^ K (buMill + ^a„7Vi?l) , (A14) 

n=l 

where a„ and 6„ are coefficients to be determined from the boundary conditions. The total 
electromagnetic wave in region 4 is therefore given by the sum of the incoming and the 
scattered fields. 

The electromagnetic fields in the other three regions are given by: 
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Using the boundary conditions: 



Ej+i X Cr, 

— * ^ 



(A18) 
(A19) 



at each interface, where indicates the radial unit vector which is perpendicular to the 
boundary, we obtain a system of 12 coupled equations to be solved. This system of equations 
can be split in two sets of six equations that we write in matrix form Ai ■ Xi — hi with for 
the first set of equations 



A, 






hn^ {kic) 









in (he) 
djn (he) 
jn ihb) 






Vn (he) 

dyn (he) 
Vn ihb) 



^dynihb) 



jn iha) 
djn {k2a) 





-jn ihb) 



Vn {ha) 
dyn (ha) 





-Vn {hb) 



-djnihb) -dyn{k2b) 



-jn {ha) 
■^.dj^iha) 






(A20) 
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and 



bi 






jn (he) 

^djn ihc) 








-.(1) 

,(2) 

L'n 

Xl - (1) 

.9n] 



(A21) 



And for the second set of equations 



A, 













djn iha) dyn {k2a) 

jn (ha) Vn (ha) 



-djnihb) -dyn{k2b) 

I&3n{hb) l^^Vnihb) -J^ihb) -yn{hb) 






djn {he) 

jn ihc) 

djn ihb) 






dyn ihc) 

Vn (he) 

dyn ihb) 



-djn (ha) 













bn 







Ujn 


djn (he) 


, X2 = 


Ujn 




M) 

Jn 







f(2) 
Jn 







. ^n _ 



(A22) 



(A23) 



These equations are solved by matrix inversion. Hence the total solution in every region is 
known. The optical cross section can now be calculated by: 



Qscatt 
Qabs 
Qtot 



2tt 



^ n=l 

27r °° 



5^(2n + l) (|a„|' + |6„|') 



(2n + 1) (|a„|^ + \brS^ - Re (a„ + 



n=l 



27r 



^ (2n + 1) Re (a„ + bn) 



(A24) 
(A25) 
(A26) 



n=l 



where a„ and 6„ are the coefficients of the scattered wave as introduced in equations (lAlSp 
and f |y414p . The wave vector k can be calculated from 



; 2 2 



(A27) 



with Ei de permittivity and /ij the permeability of the i^^"^ region. 
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